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Topological changes of the Fermi surface under pressure may cause anomalies in the low- 
temperature elastic properties. Our density functional calculations for elemental Osmium evidence 
that this metal undergoes three such Lifshitz transitions in the pressure range between 70 GPa 
and 130 GPa. The related elastic anomalies are, however, invisibly weak. The critical pressures 
considerably exceed the values for recently measured and calculated anomalies in the pressure (P) 
dependence of the hexagonal c/a lattice parameter ratio close to 25 GPa. We demonstrate that the 
latter anomalies are statistically not significant and that (c/a)(P) can be fitted equally well by a 
smooth dependence. 

PACS numbers: 71.18.+y, 71.15.Nc, 64.30.+t 



I. INTRODUCTION 

Since Lifshitz' seminal publication, 1 the search for elec- 
tronic topological transitions (ETT) of the Fermi surface 
(FS) has been a subject of permanent interest, for re- 
views see Refs. 2 and 3, As a true phase transition it 
may appear only at T = in highly pure samples, other- 
wise it is just a crossover. While alloying may change the 
band filling and, thus, the number or connectivity of the 
FS sheets, it unavoidably smears the electron states and 
the observable effects. Magnetic field can cause ETT as 
well, 4 but the field strength required to achieve a Lifshitz 
transition in a conventional metal is beyond the present 
technical possibilities. Thus, low-temperature pressure 
experiments are most promising to find measurable ef- 
fects related to ETT^ 

All thermodynamic and transport properties of a metal 
are influenced by the topology of the FS to some ex- 
tend. In the past, transport properties were in the focus 
of research, since the impact of ETT on the transport is 
stronger than on, e.g., elastic properties^ As an example, 
Godwal et alJl found a Lifshitz transition in Auln2 , where 
measurements yield anomalies in the electrical resistivity 
and in the thermoelectric power in the 2 — 4 GPa pres- 
sure range, but angle-dispersive x-ray-diffraction mea- 
surements do not indicate any structural anomaly. Very 
recently, an elastic anomaly was observed in YC05 under 
pressure, caused by an exceptionally strong van Hove sin- 
gularity in the electronic density of states (DOS)£ 

The question whether Lifshitz transitions are visible 
in the c/a ratio of the hexagonal lattice parameters of 
Zinc or if their effect is below the limit of detectability is 
controversially discussed in the literature. Earlier exper- 
imental evidence^ for an anomaly was later traced back 
to non- hydrostatic pressure conditions^ and earlier the- 
oretical confirmations of the elastic anomaly were found 
to be caused by insufficient k-point sampling^ This dis- 
cussion has recently been resumed 10 by an alternative 
presentation of the hydrostatic pressure data from Ref. 

A similar puzzling situation is present in the case of the 
5d element Osmium, which is one of the densest elements 



and the metal with the largest bulk modulus. Occelli et 
al measured the lattice parameters of Os under pressure, 
found a discontinuity in the first pressure derivative of the 
c/a ratio at about 25 GPa, and assigned this anomaly 
to a Lifshitz transition^ On the other hand, Takemura 
did not infer an anomaly from his related data obtained 
by similar measurements Subsequently, two theoreti- 
cal papers were published confirming the existence of an 
anomaly in c/a versus pressure at about 10 GPa 1 ^ and 
at 27 GPa, 14 respectively. However, no reason for this 
anomaly was disclosed in the electronic structure. In 
particular, the FS topology was found unchanged up to 
a pressure of 80 GPa* 1 ^ 

As an attempt to get a more detailed understanding of 
the dependence of the c/a ratio of Osmium under pres- 
sure, we have carried out very precise density functional 
calculations applying the FPLO code^ considering pres- 
sures up to 180 GPa. We identified three ETT in the high 
pressure region. Further, we confirm the finding of Ma 
and co-workers 1 ^ that there is no ETT within the pa- 
rameter range they considered. Our calculations yield 
pressure dependent lattice parameters that are in accord 
with those presented by the previous authors. Statisti- 
cal analysis of both our and previously published data 
eventually shows that there is no reason to assume a dis- 
cernible anomaly in the dependence of c/a on pressure 
from ETT. 



II. DETAILS OF CALCULATION 

Calculations were performed for the hexagonal close- 
packed structure (hep, P63/mmc, space group 194) which 
is the equilibrium structure of Osmium in the considered 
pressure ranged The four-component relativistic version 
of the full potential local orbital (FPLO) band structure 
code^ release 5, was used. It accounts for kinematic rel- 
ativistic effects including spin-orbit coupling to all orders. 



In the case of Osmium, the d 3 / 2 - ^5/2 splitting is 1 eV. 
Therefore it cannot be neglected in the present analysis. 
After test calculations with different basis sets we de- 
cided to use a minimum basis consisting of Afbsbp states 
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in the semicore and 6s6p5d states in the valence. All 
lower lying states were treated as core states. A k-mesh 
subdivision of 48 x 48 x 48 k-points in the full Brillouin 
zone (BZ) was used, which corresponds to 5425 k-points 
in the irreducible part of BZ. This mesh is significantly 
finer than the meshes used by Ma et al, 16 x 16 x 12 
in the BZ^ and by Sahu and Kleinman, 84 points^ in 
the irreducible part of BZ. For calculating the density of 
states, the k-mesh was refined to 96 x 96 x 96 k-points 
in the full BZ (40033 k-points in the irreducible BZ). 
The k-space integrations were carried out with the lin- 
ear tetrahedron methodJ^ For the exchange-correlation 
potential we employed the local density approximation 
(LDA) in the version proposed by Perdew and Wang, 
1992^ 



III. TOTAL ENERGY CALCULATIONS 

The upper panel of Figure [T] shows the total energy 
of Osmium plotted against volume V, where the c/a 
ratio has been relaxed for each volume, E tota i(V) = 
mhi( c / a ) E tota \(V, c/a), which is the correct condition for 
hydrostatic pressure. The red (gray) dashed line shows 
the experimental volume (taken from Ref. [ill ), deter- 
mined in a room-temperature experiment. The calcu- 
lated equilibrium volume is 0.3% smaller than the mea- 
sured room-temperature volume. Accounting for thermal 
expansion provides an almost perfect agreement. It is a 
known but not yet understood feature of LDA, that it 
reproduces the ground state volumes of the heavy 5d el- 
ements much better than gradient approximations^ We 
take the agreement in the present case as justification to 
rely on LDA in the further investigations. 

The inset of Figure [TJ upper panel, proves the excep- 
tional stability of the calculations. Each data point in 
the figure is taken from independent self-consistency runs 
and independent c/a optimizations. The FPLO codei^ 
yields a very smooth energy-curve, with numerical noise 
well below the /iHartree range. We took advantage of 
this stability in the evaluation of the equation of state 
(EOS) and the pressure dependence of c/a, presented in 
Section EH While the calculations yield directly the two 
dependences E to tai(V) and (c/a)(V), in experiment the 
pressure P = —(dE tota i/dV) may be measured indepen- 
dently. Instead of relying on an analytic fit of Birch- 
Murnaghan type for E tota i(V) we calculated the deriva- 
tive numerically by means of the three-point formula, 
f'((xi + x 2 )/2) ~ (f(x 2 ) - f(xi))/ (x 2 - xi), considering 
neighboring V values. This seemingly less sophisticated 
procedure solves a dilemma, which is present for analytic 
fits. If an analytic fit is used for the whole volume range 
considered, the anomaly will be removed by the fit. If the 
analytic fit is done piece- wise, anomalies will necessarily 
be produced at the end points of the respective fit regions, 
due to different noise in the adjacent regions. Previous 
authors used one fit for the whole range i 13 i 14 We checked 
our procedure by comparing the evaluated pressures with 
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FIG. 1: (color online). Upper panel: total energy versus unit 
cell volume with relaxed c/a ratio. The red (gray) dashed 
line marks the experimental room-temperature volume. Inset: 
magnification of the total energy versus volume curve around 
the minimum. In the large graph this region is marked with a 
green (gray) box. Lower panel: equation of state of Osmium. 
Blue (black) circles: present calculations; green (gray) pluses 
and red (gray) crosses: experimental data from Occelli et a/.— 
and from Takemura,— respectively; black line: calculations by 
Ma et a/.—. The green (gray) dashed-dotted line denotes the 
volume where Occelli et a/.— find an elastic anomaly. The 
blue (black) dashed lines indicate the volumes where Lifshitz 
transitions are found in the present calculations. 



five different fits of Birch-Murnaghan type. For 80% of 
the pressure range between zero and 180 GPa, our nu- 
merical fit lies within the P(V) values spanned by the 
different analytic fits. The maximum difference between 
any of the analytic fits and the numerical fit amounts to 
0.2 GPa, and the mean difference amounts to 0.1 GPa. 

Figure [TJ lower panel, shows the EOS of Osmium. 
Green (gray) pluses and red (gray) crosses denote ex- 
perimental data by Occelli et ali^ and by Takemura^ 
respectively. Both data sets do not deviate from each 
other on the scale of this figure. The blue (black) cir- 
cles denote results of the present calculations, obtained 
from the E(V) data, upper panel of Figure [TJ by nu- 
merical differentiation. The EOS calculated using LDA 
coincides with the experimental data in the low-pressure 
range up to about 40 GPa. It slightly over-estimates the 
pressure in the higher-pressure range. 
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Results from calculations by Ma et ali^ are depicted 
by a black line. These calculations were carried out in 
the GGA which systematically over-estimates the zero- 
pressure volume of heavy 5d metals^ Accordingly, the 
GGA data run almost parallel to the LDA data, with a 
volume offset of about 2%. An alternative presentation, 
P(V/Vo) with Vb = V{P = 0), yields a very good coinci- 
dence of both calculated data sets. As stated before, the 
offset explains the fact that Ma et al. did not find any 
ETT in their calculations up to 80 GPa^ while we find 
the first ETT at about 72 GPa. 

In the lower panel of Fig. [TJ there is no anomaly of the 
EOS visible in any data set at any pressure, in agreement 
with the earlier presentations of the quoted data and with 
the weakness of the Lifshitz transitions disclosed in this 
paper. 



IV. THE FERMI SURFACE UNDER PRESSURE 

It was suggested by Occelli et al, that the observed 
anomaly in (c/a)(P) should arise from an ETT. Thus, we 
first investigate the volume dependence of the FS in the 
range 22.0 A 3 < V < 28.0 A 3 . This volume dependence 
translates into a pressure dependence through the EOS 
given in the lower panel of Figure HJ P < 180 GPa. As in 
the evaluation of the total energy, c/a has been relaxed 
for each volume. 

At zero pressure, the FS of Osmium has four sheets, 
see Figure [2j The first sheet consists of small hole el- 
lipsoids located at the line L-M. Occelli et al. suspected 
these ellipsoids to disappear under pressure.— The sec- 
ond sheet is a big monster surface, which is closed in 
z-direction (compare Figure [3]) but connected in the x-y 
plane. The third and fourth sheets are closed surfaces, 
which are nested and centered around the T-point. The 
inner one of the two is waisted. Our calculated zero- 
pressure FS agrees well with measurements of Kamm 
and Anderson^ with the calculated FS of Smelyansky 
et al., 21 and with the more recent calculation of Ma et 

Figure [3] shows a comparison of the FS at zero pressure 
(first row, V = 27.87 A 3 , c/a = 1.5850) with the FS cal- 
culated at a very high pressure of about 180 GPa (second 
row, V = 22.0 A 3 , c/a = 1.5968). It is obvious, that the 
small hole ellipsoid between L and M does not disappear 
under pressure. On the contrary, it grows. Neverthe- 
less, we find a Lifshitz transition in the first FS sheet: 
a further tiny hole ellipsoid appears under pressure at 
the L-point. Another Lifshitz transition takes place in 
the second FS. Here, a neck is created under pressure, 
centered at L. The third FS sheet is not displayed in Fig- 
ure [3] because its topology is preserved in the considered 
pressure range. Finally, a third ETT is found inside the 
fourth FS sheet, where a hole ellipsoid appears under 
pressure at the T-point. 

At which pressure/ volume the FS changes its topology 
depends on the c/a ratio. This can be seen in Figured 




FIG. 2: (color online). The Fermi surface of Osmium at zero 
pressure. First, second, third, and fourth sheets are given 
separately in the upper left, upper right, lower left, and lower 
right panel, respectively. Red (black) and yellow (light-gray) 
colors denote occupied and unoccupied sides of the Fermi sur- 
face, respectively. 



Reducing the volume under hydrostatic pressure, start- 
ing at the zero-pressure value, one finds the following 
sequence of transitions: At first, the additional ellipsoid 
at the T-point appears at V = 24.60 A 3 , then the neck 
is formed at V = 24.20 A 3 , and finally the additional 
ellipsoid at the L-point appears at V = 23.20 A 3 . The 
related transition pressures are 72 GPa, 81 GPa, and 122 
GPa. Worth mentioning, a peculiarity relates the latter 
two ETT. Namely, two necks are formed at first at the 
mentioned volume which are situated at the line L-H. 
These necks merge at L at the same pressure that lets 
the L-point ellipsoid appear. This coincidence is due to 
a degeneracy of the first and second sheet bands along 
L-A, see Figure [5] 

If we presume the experimentally measured and ex- 
trapolated (c/a)(V), the ETT take place at pressures of 
about 50 GPa (ellipsoid at the T-point), at about 87 GPa 
(neck at the line L-H) and at about 122 GPa (ellipsoid 
at the L-point). 

The discussed topological changes can be clearly iden- 
tified in the band structure plots. For each hole ellipsoid 
appearing under pressure one expects a related maximum 
in the band dispersion to cross the Fermi energy. In the 
case of the neck, a saddle-point of band dispersion should 
cross the Fermi level. This saddle-point does not lie on a 
symmetry-point but only on a symmetry-line of the Bril- 
louin zone and thus appears only as a maximum in the 
band structure plot. 

Figure [5] shows the band structure for zero pressure 
(upper panel) and under a high hydrostatic pressure of 
about 180 GPa (lower panel). The loci of the three dis- 
cussed ETT are marked by red (gray) full-line circles. 
The pressure-induced Fermi level crossing of all these 
three band maxima is obvious. The blue (black) dashed- 



line circles in Figure [5] mark the place of the disappear- 
ance of a hole pocket proposed in Ref. [ll|. The related 
maximum does, however, not fall below the Fermi level 
under pressure but gets shifted to higher energy instead. 

Summarizing this section, we find three distinct Lif- 
shitz transitions in the pressure range between 70 GPa 
and 130 GPa. In the range of 25 GPa the present cal- 
culations do not yield an anomaly of c/a caused by the 
electronic structure. None of the transitions disclosed 
here were found in the earlier calculations restricted to a 
pressure range just below the first observed transition^ 
Note, that the Fermi surface topology in the latter pub- 
lication was investigated up to 80 GPa. The related 
volume, however, was slightly larger than our predicted 
first critical volume. The reason is that Ma et al. used 
the generalized gradient approximation (GGA), while we 
used LDA. 



V. DENSITY OF STATES 

ETT should manifest themselves in van Hove singular- 
ities of the DOS passing through the Fermi level under 
pressure. The upper panel of Figure [6] presents the DOS 
of Osmium at zero pressure and at about 180 GPa. The 
5d band is considerably broadened at such a high pres- 
sure, in comparison with the zero pressure case. How- 
ever, the Fermi level is situated in a smooth valley of 
the DOS in the whole pressure range. In particular, no 
Fermi level crossing of any singularity is obvious in this 
overview panel. 

The lower part of Figure [6] shows the DOS close to the 
Fermi level in about hundredfold magnification, for five 
different pressures including the three transition pres- 
sures and the two limiting cases of the upper panel. A 
tiny anomaly, caused by the Lifshitz transitions, is found 
moving through the Fermi level. Despite the huge num- 
ber of 96 3 k-points in the full BZ used for these calcula- 
tions, the estimated DOS resolution that can be achieved 
for the present band dispersions by means of the linear 
tetrahedron method amounts to 50 meV which is about 
the width of the observed anomaly. Thus, no details of 
this feature can be resolved, it remains unclear to which 
extent the individual ETT contribute to it. (Note, that 
the shape is changing with pressure. This fact also hints 
to resolution problems.) 

As the ETT anomaly amounts to about 10 -3 of the 
total DOS, it is probably impossible to detect it in the 
elastic properties. This becomes even more clear if we 
consider the volume dependence of the DOS at the Fermi 
energy, Figure [7J which decreases almost linearly with 
unit cell volume. The red (gray) solid line shows a linear 
fit to the data. No peculiarity is visible at the Lifshitz 
transitions, marked by blue (black) dashed lines. To un- 
derstand this point, we have estimated the strength of 
the van Hove singularity related to the first ETT, the 
appearance of the ellipsoid at T. The extra contribution 



to the DOS due to this ellipsoid can be calculated as 

ir / x. V ISm x m y m z r 

= ^2 Y ff V(^cr-£) (1) 

with V = 27.87 A 3 . Here, h 2 /2mi denotes the 
curvature of the ellipsoid in direction z, which has 
been fitted to h 2 /2m x = 171.4 eVA 2 , h 2 /2m y = 
161.4 eVA 2 and H 2 /2m z = 97.1 eVA 2 . Then Sg(e) « 
0.0009 eV _ V( £ cr - e)/eV or about 0.00009 eV" 1 at 10 
meV distance from the position of the van Hove singu- 
larity, e cr . Here, the crucial point is the tiny prefactor 
of the square-root singularity. This extra contribution to 
the pressure of the Fermi gas will not be visible in elas- 
tic data obtained by today's experimental or numerical 
means. 

Summarizing this section, we find that the disclosed 
Lifshitz transitions are very weak and probably not de- 
tectable by measurements of elastic properties. How- 
ever, there might be a chance of observing anomalies in 
magneto-transport or thermopower under pressure. 

VI. c/a RATIO UNDER PRESSURE 

We now proceed to a comparison of several measured 
and calculated variants of the pressure and volume de- 
pendence of the c/a ratio. A critical evaluation and data 
analysis will help to clarify the question, if there is any 
anomaly in the data close to 25 GPa (or 10 GPa), as 
proposed by several authors j n i 13 i 14 

Figure [8] shows the dependence of c/a on volume and 
on pressure. Similar to the EOS case, both experimen- 
tal data sets coincide within the error estimates given 
in Refs. [ll| and [H. The three computed data sets run 
more or less parallel with the experimental ones. Un- 
like the EOS data, the offset of the LDA results of about 
0.005 (blue (black) circles, present calculation, and vio- 
let (black) triangles, calculation by Sahu and Kleinman) 
is larger than the offset of the GGA results by Ma et 
al. (black diamonds). Additionally, we have introduced 
another data set (blue (black) open circles) in the lower 
panel of Figure [HI denoting our calculated data shifted 
down by 0.00515. These offset-corrected LDA data show 
a remarkable coincidence with the measured data both 
in slope and curvature. Thus, they can be used to ex- 
trapolate the experimental data into the high-pressure 
region in order to estimate experimental transition pres- 
sures given in Section HVl 

In (c/a)(V/Vo) no anomaly is visible at any of the Lif- 
shitz transition volumes (blue (black) dashed lines), and 
no anomaly is visible either at a volume of about 0.95^o, 
corresponding to 25 GPa. The proposed anomalies in 
(c/a)(P) will be analyzed subsequently. Before, we check 
the dependence of c/a on the volume in detail, which 
looks linear at first glance. 

We performed linear fits to the individual data sets 
in the volume range 0.87 < V/Vq < 1, accounting for 
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Xn 


a 


H/A 3 


experiment (Occelli et a/.) 1 ^ 


0.73 


7.4 


10~ 5 


27.949 


experiment (Takemura)^ 


1.22 


1.2 


io- 4 


27.977 


present calculation 




5.9 


10~ 5 


27.87 


calculation by Ma et al— 




8.1 


10~ 5 


28.48 



TABLE I: Quality of linear fits to the (c/a)(V/V ) data dis- 
played in Figure [U upper panel, in the range 0.87 < V/Vq < 1. 
For the definition of Xn and a see text. The fourth column 
displays the values used for Vo- Values of adjustable parame- 
ters are given in the Appendix. 

the known uncertainties of the experimental values, and 
considering the theoretical data points with equal weight. 
To char acterize the q uality of the fits, Table [J displays 
Xn = \/x 2 / ( n p — n e) where x 2 denotes the sum of the 
error weighted square deviations, n p the number of the 
data points, and n e the number of adjustable parameters, 
here 2. Moreover, for all data sets, TableUgives estimates 
of the standard deviation a presuming validity of the 
linear dependence and uniform uncertainties of all data 
points included. Two conclusions can be drawn from 
this check: (i) the dependence of c/a on volume is linear 
within experimental accuracy according to the Xn values, 
and (ii) the random errors of the theoretical data are 
comparable with the experimental inaccuracy. A linear 
fit of our calculation results for the whole volume range, 
0.78 < V/V < 1, yields a a value of 8.7 • 10~ 5 , resulting 
from a slight non-linearity at higher compressions. 

Now we come to the key question of a possible anomaly 
in (c/a)(P) at about 25 GPa, green (gray) dashed-dotted 
line in the lower panel of Figure [HJ We stress again, in 
accordance with Refs. [ll| and [3, that there is no doubt 
about the linearity of (c/a)(V/Vo) within experimental 
precision. Thus, if (c/a)(P) would be piecewise linear, 
as suggested by Occelli et al., Sahu and Kleinman, as 
well as Ma et ah, then P(V/Vo) had to be piecewise lin- 
ear as well. In turn, this would imply a zero pressure 
derivative of the bulk modulus, B' = 0, in each of the 
pressure regions below and above the anomaly - a clear 
contradiction with the EOS fits by Occelli et at. who find 
B' « 4 in both regions*^ 

Given a non-zero pressure derivative of the bulk mod- 
ulus, P(V) could most simply be approximated by a 
second-order polynomial. Since c/a is linear in the vol- 
ume, this approximation translates into a square-root de- 
pendence of c/a on pressure: 

(c/a) (P) = (c/a) m + dy/(P + P m ) , (2) 

where (c/a) m , d, and P m are fit parameters. 

It is clear, that a square-root behavior can be fitted 
by two linear pieces within a restricted range. Table [Til 
gives a comparison of square-root and piece-wise linear 
fits between zero and 60 GPa, which is the range covered 
by all of the data sets. (We do not take into account the 
data by Sahu and Kleinman, since they have only four 
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piece-wise linear function 








kink: 25 GPa 


kink: 27 GPa 




Xn 
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Xn 
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Xn 


a 


exp. (Occelli et al.)^ 


0.69 


6.8 • 10" 5 


0.71 


7.4 • 10" 5 


0.85 


8.5 • 10" 5 


exp. (Takemura) 1 ^ 


1.27 


1.3 • 10" 4 


1.28 


1.3 • 10" 4 


1.26 


1.3 • 10" 4 


present calculation 




7.5 • 10~ 5 




1.1 • 10" 4 




1.0 • 10" 4 


calc. by Ma et al^ 




7.8 • 10" 5 




1.0 • 10" 4 




9.7- 10~ 5 



TABLE II: Quality of different fits to the data displayed in 
Figure [HI lower panel, in the pressure range between zero and 
60 GPa. For the definition of Xn and a see text. Values of 
adjustable parameters are given in the Appendix. 



points in the considered pressure range.) According to 
Refs. [ll| and [3 we assumed the kink between the linear 
pieces at the suggested 25 GPa and 27 GPa, respectively. 

Table [II] shows that both kinds of fits work almost 
equally well, with a small advantage of the square-root 
fit. Additionally for the data by Occelli et a/., we per- 
formed a four parameter piece- wise linear fit, allowing 
the critical pressure P c to vary: Only a slightly better 
approximation was obtained (% n = 0.51 for the complete 
data set), but P c amounts to 19 ±2 GPa. Finally we have 
checked parabolic and cubic fits, which yield very similar 
Xn and <j values compared to the square-root fit. 

According to Tab. [Til and the non- vanishing pressure 
derivative of the compressibility argument above, we be- 
lieve that there is not enough evidence to conclude an 
anomaly at about 25 GPa. Regarding the DFT calcula- 
tions, no peculiarity in the band structure or Fermi sur- 
face topology is observed in this pressure range. Thus, 
it would be artificial to expect an anomaly in the com- 
puted elastic data. Of course, regarding the experiments, 
other than electronic origins of an anomaly cannot be ex- 
cluded. However, in our opinion, such an anomaly is not 
proved by in the available data sets. Error correlations 
mentioned in Ref. [ll| could be caused by the change of 
the experimental setup at 26 GPa. 



VII. SUMMARY 

We predict three Lifshitz transitions to occur in hexag- 
onal close-packed Osmium under pressure up to 180 GPa. 
The corresponding van Hove singularities in the density 
of states are probably too small to produce any measur- 
able effect in the elastic properties but may be detectable 
in transport properties. Since the lowest critical pres- 
sure amounts to 70 ... 80 GPa, it is not related with a 
previously suggested kink at 25 GPa in the pressure de- 
pendence of c/a. We demonstrate that this kink can be 
replaced by a smooth dependence without any change in 
the statistic significance. Thus, we cannot confirm the 
electronic origin of the proposed anomaly. 
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A 


10 4 5/GPa- 1 


lO^/GPa- 1 


exp. (Occelli et a/.)ii 
exp. (Takemura)^ 


1.579933(9) 
1.57995(4) 


1.23(2) 
1.17(3) 


0.80(3) 
0.86(3) 


present calculation 
calc. by Ma et ali^ 


1.58525(2) 
1.57729(7) 


1.117(8) 
1.46(3) 


0.955(6) 
1.06(4) 



APPENDIX 



1. Linear fit to (c/a)(V/V ) 
c/a = A- (V/Vo) + B 





A 


B 


exp. (Occelli et al.)^ 
exp. (Takemura)^ 


-0.0539(4) 
-0.0540(6) 


1.6339(4) 
1.6339(6) 


present calculation 
calc. by Ma et ali^ 


-0.0594(2) 
-0.0678(6) 


1.6445(1) 
1.6450(6) 



2. Square-root fit to (c/a)(P) 





(c/a) m 


d/GPa- 1 / 2 


Pm/GPa 


exp. (Occellie£ al.)^ 
exp. (Takemura)^ 


1.572(1) 
1.572(2) 


0.00146(9) 
0.0015(1) 


28(5) 
30(9) 


present calculation 
calc. by Ma et ali^ 


1.5753(8) 
1.5708(10) 


0.00162(5) 
0.00157(8) 


37(4) 
16(4) 



Kink at 27 GPa 



A - c/a at P = OGPa 

B - slope between GPa and 27 GPa 

C - slope between 27 GPa and 60 GPa 





A 


lO^/GPa" 1 


10 4 C/GPa" 1 


exp. (Occelliet aZ.)il 
exp. (Takemura)^ 


1.579935(2) 
1.57995(4) 


1.219(14) 
1.16(3) 


0.78(4) 
0.84(3) 


present calculation 
calc. by Ma et ali^ 


1.58525(2) 
1.57730(6) 


1.113(8) 
1.45(3) 


0.940(7) 
1.02(4) 



3. Piece-wise linear fit to (c/a)(P) 
a. Kink at 25 GPa 



A - c/a at P = OGPa 

B - slope between GPa and 25 GPa 
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FIG. 3: (color online). Comparison between the Fermi surface 
of Osmium at zero pressure (V — 27.87 A 3 first row) and un- 
der high pressure (second row, P « 180 GPa, V = 22.00 A 3 ). 
The upper left panels show the first FS sheet, the upper right 
and the lower left panels show the second FS sheet in a side 
view and in a top view, respectively. In order to make the con- 
nectivity better visible, the side view is shifted in ^-direction, 
such that the BZ edge lies in the middle of the picture. The 
lower right panels display a cut through the fourth sheet of 
the Fermi surface in the y-z plane. Colors are explained in 
Figure [2] 
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FIG. 4: (color online). Dependence of the FS topology on 
volume and c/a. Black symbols with error bars: calculated 
relaxed c/a versus unit cell volume. Green (gray) full line: 
ETT at the T-point. Red (gray) dashed line: ETT at the 
symmetry-line LH. Blue (black) dash-dotted line: ETT at 
the L-point. 



Os, V=27.87A 3 , c/a=1.5850 




Os, V=22.00A 3 ,c/a=1.5968 
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FIG. 5: (color online). Band structure of Osmium. Upper 
panel: zero pressure; lower panel: P w 180 GPa. The red 
(gray) circles (full lines) mark points in k-space where a band 
maximum crosses the Fermi level under pressure. The blue 
(black) circle (dashed lines) marks the location of an ETT 
proposed by Occelli et a/..— 
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FIG. 6: (color online). Upper panel: Total density of states of 
Osmium at two different volumes. Maroon (black) full lines: 
zero pressure; red (light-gray) dashed lines: pressure of about 
180 GPa. The Fermi level is situated in a smooth valley of 
the DOS. Five lower panels, from left to right: magnification 
of the DOS around the Fermi level at zero pressure; at that 
pressure, where the ellipsoid at the T-point appears, at that 
pressure, where the neck is created, at that pressure, where 
the ellipsoid at the L-point appears, and at a pressure of about 
180 GPa. 
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FIG. 7: (color online). DOS at the Fermi energy versus vol- 
ume. The red (gray) solid line shows a linear function fitted to 
the data. The blue (black) dashed lines mark those volumes, 
where the Lifshitz transitions take place. 
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FIG. 8: (color online). Upper panel: relaxed c/a-ratio against 
normalized volume; lower panel: relaxed c/a-ratio against 
pressure. The meaning of the symbols +, x, • and green 
(gray) dashed-dotted line is the same as in Figure [1] lower 
panel. Moreover, black diamonds denote calculations by 
Ma et a/..— Additional data from calculations by Sahu and 
Kleinman are denoted by violet (black) triangles in the lower 
panel.— The blue (black) open circles are a shift of our data 
down by 0.00515. Error-bars are below symbol size. The 
position of the Lifshitz transitions (blue (black) lines) in the 
lower graph was obtained from the FPLO EOS. 



